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We consider stochastic energy balance and entropy production (EP) in a generalized Langevin 
dynamics of macrospins, allowing for both amplitude and direction fluctuations, under external 
magnetic field. EP is calculated using Fokker-Planck equation, distinguishing between reversible 
and irreversible parts of probability currents. The system entropy increases due to irreversible 
non-equilibrium processes, and reduces as heat dissipates to surrounding environment. Using path 
probability distributions of time-forward trajectories and conjugate trajectories under time reversal, 
we obtain fluctuation theorems (FT) for total stochastic EP. We show that the choice of conjugate 
trajectories is crucial in obtaining entropy like quantities that obey FTs. 

PACS numbers: 05.40.-a, 05.40.Jc, 05.70.-a 


I. INTRODUCTION 


Stochastic spin dynamics under magnetic fields and 
the influence of other spins, plays an important role in 
understanding magnetic properties of condensed matter 
systems. With the advent in nano-technology, the size of 
magnetic devices like magnetic read head and random ac¬ 
cess memory are being reduced consistently This makes 
them vulnerable to thermal fluctuations [l|, [^. Under¬ 
standing the role of stochasticity in such devices is thus 
becoming important, even from the perspective of better 
control of their performance The classical dynam¬ 

ics of a magnetization m under external field H is de¬ 
scribed by the Heisenberg motion rh = ym x H where 
7 denotes the gyromagnetic ratio This dynamics, 
evidently, conserves the amplitude m = |m|. When cou¬ 
pled to a heat bath, the dynamics gets stochastic and 
is often e:™ressed as a stochastic Landau-Lifshitz (sLL) 
equation [91 

rh = ym X [(H + h(t)) — r/{m x H)]. (1) 

Kubo and Hashitsume argued for the introduction of dis¬ 
sipation term —rj'm x (m x H) along with the stochas¬ 
tic fluctuation h(t) to obtain fluctuation-dissipation rela¬ 
tion 0. Here h(t) is regarded as a Gaussian white noise 
with (h(f)) = 0, and (h(t)®h(f')) = 2Z?gli5(t — t') with 1 
denoting an identity matrix, = rj'ksT/V where T is 

the temperature, V the physical volume of the macrospin, 
ks Boltzmann constant. The sLL equation was inde¬ 
pendently derived using the Zwanzig formalism of cou¬ 
pling the spin dynamics with harmonic bath and tak¬ 
ing the Markovian limit M- This equation describes a 
stochastic rotational dynamics of magnetization, keeping 
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the magnitude m conserved. In Ref. ll|, the proper¬ 
ties of Fokker-Planck equation corresponding to a related 
Landau-Lifshitz-Gilbert (LEG) equation were analyzed 
in detail. The constant m dynamics is a good approxima¬ 
tion for bulk ferromagnets at room temperature, where 
the transition temperatures for ferromagnetic to para¬ 
magnetic phase transitions are much larger. Stochastic 
fluctuations of m occur within a ferromagnetic domain 
due to exchange interactiom with enhanced effect near 
the transition temperature [8]. The fluctuations in m be¬ 
comes dominant in bulk ferromagnets only at high tem¬ 
peratures. On the other hand, due to enhanced relative 
fluctuations in small ferromagnetic domains, e.g., in a 
macrospin, the transition temperatures get largely sup¬ 
pressed with reduction of system size [l^, 13 1, enhancing 
the fluctuations in m even at room temperature. Re¬ 
cently, a generalized Langevin spin dynamics has been 
proposed that captures longitudinal fluctuations in the 
spin magnitude, as well as the stochastic rotation dy¬ 
namics of its orientation [l^ . In this paper, we present 
stochastic thermodynamics of a macrospin system, de¬ 
riving stochastic energy balance relation and fluctuation 
theorems for probability of entropy production. 


During the last two decades, a theoretical descrip¬ 
tion of stochastic thermodynamics has been developed 
to describe non-equilibrium small systems having en¬ 
hanced relative fluctuations, using stochastic counter¬ 
parts of thermodynamic variables like energy, work, en¬ 
tropy etc. [isl - l^ . While the possibility of second 
law violating stochastic trajectories was recognized long 
back [ 2 ^, it took several decades before it was shown 
that probabilities of such trajectories in steady state 
P{—Ast), with —Ast denoting negative entropy pro¬ 
duction, are exponentially suppressed with respect to 
the positive entropy producing ones via the relation 
P{Ast)/P{—Ast) = exp{Ast/kB) [23 - 1^ . This relation 
is known as the detailed fluctuation theorem. This, and a 
related integral fluctuation theorem (exp(—Ast/Zcs)) = 
1, which is equivalent to the Jarzynski equality for non¬ 
equilibrium transformations from an initial equilibrium 





state to a final state that eventually reaches equilib¬ 
rium, have been derived l27H^ . These theorems 

were verified using experiments on colloids [sil - l^ , gran¬ 
ular matter [sj, Isal. a nd used to obtain the free energy 
landscape of RNA |36l a and torque produced by Fl- 
ATPase motor proteins m. 

In the following section, we present the generalized 
Langevin dynamics of macrospins (GLDM). We discuss 
its motivation, and corresponding projected dynamics 
along the longitudinal and transverse directions. Next, 
we study its stochastic thermodynamics, first deriving 
the stochastic energy balance, and then entropy pro¬ 
duction using Fokker-Planck equation and ratio of time- 
forward and conjugate path-probabilities. Our analysis 
shows that it is possible to obtain fluctuation theorems 
for entropy like quantities, each of which emerges out of a 
specific way of choosing conjugate trajectories. The time- 
reversed trajectories give fluctuation theorems in terms 
of EP in reservoir given by the dissipated heat, which 
is consistent with the results of Fokker-Planck equation. 
Another possible choice of conjugate trajectories leads 
to an entropy like quantity that also involves gyroscopic 
work done due to magnetic field induced spin torque. 
This quantity also obeys both detailed and integral fluc¬ 
tuation theorems. We present discussions interpreting 
our results. Finally, we conclude by presenting a sum¬ 
mary. 


II. MODEL 


Consider a macrospin having magnetization m, and 
volume V. The GLDM for the macrospin in presence 
of a time-dependent external magnetic field H(f) can be 
written as [l4l | 


rh = 7 


m X H(t) — r/ 


dm 


h(t) 


( 2 ) 


where rh = dm/dt. The Langevin heat bath is charac¬ 
terized by the dissipation coefficient 77 , and the Gaussian 
white noise h the components of which obey (h(t)) = 0 , 
(h(t) 0 h{t')) = 2Dol6{t — t') with Dq = ■qkBT/V. In 
the above equation m x H denotes a non-conservative 
spin-torque. The energy density is given by 


g = if L- m.H(t)) 

where, /l = (3) 

is the Landau free energy density, and the effective 
magnetic field H®® = -dg/dm = -f H(t), with 
H'"* = (a — bm^)m being the mean field contribution 
due to collective spin alignment, g can be expressed as 
g = —m.H®®. Jl denotes the Landau free en ergy density 
having two equivalent minima at m = ±y/a/brh, with 
a = ao{Tc — r) > 0 in the ferromagnetic phase, where 
Tc is the transition temperature, and T is the temper¬ 
ature of the system Q. With reduction of system size. 


Tc decreases. It was shown for three dimensional Ising 
clusters with total number of spins N, the transition tem¬ 
perature decreases with reduction in macrospin size N as 
Tc - T~(l - l/TV-^) where cj) a; 1/3, and denotes the 
transition temperature of thermodynamically large sys¬ 
tem Thus for small enough size of a macrospin, Tc 
approaches T from above, thereby increasing fluctuations 
in TO. The term —m.H in g is due to external magnetic 
field H, and shifts the global minimum towards positive 
rh. Thus the GLDM may be expressed as 

m= [m X H(t) + 77 H®®(t)+h(f)] , (4) 

absorbing 7 into the definition of time, t —>■ 'yt. 

Eq. ([2]) may be motivated by drawing parallel to 
Langevin equations of motion of driven diffusing parti¬ 
cles M- Note that such a particle in one dimension 
obeys v = f(t) — gv + where v denotes the particle 
velocity, f{t) a time-dependent external force. The vis¬ 
cous dissipation —gv and the Gaussian white noise ^ are 
forces due to coupling to the heat bath, with (^(t)) = 0 , 
and = 2gkBT5{t — t'). The viscous dissipation 

—gv may be rewritten as —gdyH given that the veloc¬ 
ity dependence of the Hamiltonian Td is v‘^/2. Using this 
as a guiding principle, one may replace v by the mag¬ 
netic moment m, as both are odd parity variables under 
time reversal. Similarly, the force / may be replaced by 
the external torque due to the magnetic field m x H. 
Replacing Hamiltonian T-l by energy density g = —m.H 
for a single spin, and dyH by dg/dm = —H we obtain 
the GLDM for a single spin rii = [m x H -|- gH. + h(t)]. 
Note that in this equation, the term gll + h denotes the 
force and torque due to the heat bath [I^. Extending 
this argument to a macrospin containing large number of 
spins, one obtains Eq.(|l]) by using g = —m.H -1- Jl- In 
Eq. dH) , the term 77 H®® -1- h denotes the force and torque 
on the macrospin due to the heat bath. Throughout this 
paper, we use Stratonovich convention while interpreting 
stochastic differential equations. 

It is possible to separate the longitudinal and trans¬ 
verse dynamics of the macrospin m. Taking longitudinal 
projection, i.e., projecting Eq. ([2]) along m = m/m we 
obtain the dynamics for the spin amplitude 

rh = [gH^\(t) + + g{am — bm^)] , (5) 

where = TO.H(t) and /i|| = m.h, with (h||(t)) = 0, 

and (/i|| (t)/i|| (f')) = 2Do6{t—t'). Glearly, in this equation, 
g{H^\ +am — bm^) + h\\ is the longitudinal force due to the 
Langevin heat bath. The corresponding Fokker-Planck 
equation is dtP(m,t) = —dmj where j = —D^dynP + 
77 [aTO — bm^ + L^||]- For a time-independent magnetic 
field, setting the dissipative current j = 0 , one obtains 
the detailed balanced equilibrium distribution, Peq{m) = 
Pq exp(—gyU/fesT) with the energy density g\\ = fh — 
H\\m. 

Subtracting the longitudinal dynamics Eq.(IS|) from 
Eq.(|3]), one obtains 

rhu = m X H(t)-I-?7Hj_(t)-f hu(t), (6) 


2 





where mj_ = m —mm, Hj_ = H —m-ff|| = —rax (mxH), 
and h_L = h — mh||. It can be shown that hj_ and hj_ = 
m X h obeys the same statistics : (hj_) = 0 = (hj_), 
(hj_(t) 0 h_L(t')) = (1 — m (g) m)5{t — t') = (h_L(t) ® 
hx(t')) M- Thus one can replace h_L by hj. in Eq.(|6|. 
The resultant equation can be expressed as 

rhx = m X (H(t) + h'(t)) — rj'm x (m x H(t)), (7) 

where h' = h/m, and r/' = rj/m^. Note that for constant 
magnitude m, rhu = rii, and Eq. 0 is then same as the 
sLL equation Eq.(IT|). The term —r]'m x (m x H) denotes 
the Gilbert damping, and m x h' denotes the stochastic 
part of the total torque imparted due to the Langevin 
heat bath. 


B. Entropy production and heat dissipation 


With P(m, t) denoting the probability of hnding a spin 
in the state m at time t, the non-equilibrium Gibbs 
entropy S = —ks f dm Pin P(m, t) suggests a defini¬ 
tion of time dependent stochastic entropy of the system 
s(t) = —kBlnP(m,t) where S = (s) denotes the ensem¬ 
ble average of stochastic entropy [^. Note that, here 
and in the rest of the paper, whenever we mention en¬ 
tropy, energy or work done, it actually means the quan¬ 
tity per unit volume. To achieve this in the following, we 
replace diffusivity of magnetization Dq hy D = riksT. 

The rate of change in stochastic entropy is given by 


s 

ks 



( 11 ) 


III. RESULTS AND DISCUSSIONS 
A. Energy conservation 

The macrospin undergoes a relaxation dynamics in 
the Langevin heat bath, settling into an average uni¬ 
directional precession around the field H. Unlike sLL 
equation which describes motion of spin constrained to 
have a constant magnitude, here amplitude m = |m| is 
not conserved, and obeys the above mentioned distribu¬ 
tion Peq(m). The rate of change of energy density g of 
the macrospin is given hy g = —m.H®® — m.H, with 
= H -b Pint and Hint = -dfL/dm. Note that H 
and Hint shares the same symmetry under time reversal, 
as does H and m. Substituting Eq. ([2]) in the ex 
of g we obtain stochastic energy balance. 


eff 


g = q + w 
where q = —m.H 
and w = —m.H . 


We used the sign convention that stochastic heat and 
work done are positive if they increase the energy of the 
system. Note that the stochastic energy balance pre¬ 
sented above, is a relation between energy density, work 
density and heat absorption per unit volume. Here q 
and w represent the rate of heat absorbed by the system 
and the rate of work done on the system respectively. 
Note that for driven diffusive particles, heat absorbed by 
the system is given by u(—ryu -|- ^), and we motivated 
the GLDM equation by replacing v by — H, and ^ by h. 
Thus it is only natural to identify —{gH^ -f h.H) as heat 
absorbed by a single spin. For a system of spins present 
in the macrospin, H has to be replaced by H®®. Thus 
q = — 4. h jj®®] = —m.H®®. This suggests a 

stochastic version of Glausius entropy production in the 
heat bath in the form —q/T. In the following, we present 
a careful analysis of entropy production. 


where the probability density P{m,t) obeys the 
Fokker-Planck equation dtP = —dm-J with dm = 
and probability flux J. Note that un¬ 
der time reversal t, m and H change sign. Thus J = 
jr®v _|_ where J''®'' is the reversible current that does 
not change sign under time reversal, and is the irre¬ 
versible current that changes sign [^. The i-th compo¬ 
nent of these currents are given by 


jr^ = 

J'r = gHfP-DdrmP. 


N,P 

reff 7 


( 12 ) 


where, the component of spin-torque due to external 
magnetic field Ni = (m x H)^. Using Eq. (fT^ in Eq. (fTTl) 
to replace ^mP^ one can express the rate of change in 
stochastic entropy as 



s 

dtP 


ks 


(8) 


dtP 

(9) 


P 

(10) 

At this stage. 

let us 




rui 


PD 
rhi 


^ zyefF • 


J] 


1 


PD 


ksT 


q- 


(13) 


(i) over trajectories, and (ii) over phase space by integrat¬ 
ing over all m with probability P(m,t). The trajectory 
average of the components of magnetization dynamics 
depends on both reversible and irreversible parts of prob¬ 
ability flux, (rhiI m, t) = Ji/P = (J[®'' -|- = Ni + 

/P [ 3 ^. Thus after the trajectory average one can re¬ 
place {Jf^mi)/{PD) by \{JrNi)/{PD) + {Jf^f/{P^D)]. 
Now let us perform averaging over the probability density 
P(m,t) by multiplying Ea. dT^ throughout by P(m,t) 
and performing integration over m. The conservation of 
probability f dmP = 1 leads to f dm dtP = 0. Thus one 
obtains the final average 


ks 


(f) 

ks 


1 

P 


dm 


[Jr? 


h 


dm jrN, - 


Ml 

keT' 


Now using the expression in Ea. dT^ one can show that 
the second term in the above equation f dm Jf^Ni = 0. 
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This term vanishes, as (i) Hf^Ni = 0 due to vec¬ 
tor identities H.(m x H) = 0 and m.(m x H) = 0, 
(ii) f dm NidruiP = 0 using integration by parts. Thus 


^ = (s) = 


VT. 


dm 




irr'i2 


P 


M 

T 


= n - 


(14) 


where II = ^ J dm ^'’'p ^ is the EP in the system due to 
irreversible processes quantified by J™, and <i> = {q)/T 
is the entropy fiux to the reservoir due to average heat 
loss. At this stage, it is interesting to note that as the 
system gets isolated from the heat bath, i.e., 77 —>• 0 , 
n ~ 77 ^ 0, a result expected for EP in an iso¬ 
lated system. The total EP in the combined system 
and reservoir obeys the second law of thermodynamics, 

St = S-\-^=:^J dm > 0, where the equal¬ 

ity denotes equilibrium with = 0. At steady state, 
<& = n = (g)/T, and average change in energy {g) = 0 
leads to {q) = — (w) = (m • H). Thus one can express the 
average entropy flux as $ = (m • H)/T. The steady state 
EP in the reservoir is due to non-equilibrium processes 
driven by time-dependent external field H(t). 

The above discussion shows that the stochastic EP in 
the reservoir is 



This quantity can be both positive or negative. Oono 
and Paniconi introduced a concept of housekeeping 
heat, which is the heat dissipated to keep the system at 
non-equilibrium steady state. As we have seen above, 
at steady state, the average heat dissipated is equal to 
the mean work done by the system (— m • H). Thus 
the expression of stochastic housekeeping heat genera¬ 
tion q/j = — m • H. This gives the rate of excess heat 
generation q^ = q —qh = m-'H. —m- H®®. If one changes 
the magnetic field from some initial value to a final value, 
the average excess heat generation remains non-zero tran¬ 
siently before the system relaxes from one steady state 
to another. 


C. Equilibrium detailed balance 

The steady state condition is given by dmAJr'' + 
JfA = 0 . At equilibrium, the dissipative current must 
vanish, = 0. This leads to the condition dP/P = 
fdHf^dmi. The equation can be integrated for a time- 
independent external field H to give 

P = Pq exp[-/3{/L - m • H}], 

where /l = —{aj2)m^ + {b/A)m^. Using the rela¬ 
tion J*’’'' = 0 in the steady state condition one obtains 
dmiJ^'' = 0, which is readily obeyed. These two rela¬ 
tions, J-''"' = 0 and drriiJr'' = 0 define the equilibrium 
detailed balance condition. A time-dependent magnetic 
field brings the system out of equilibrium, and allows EP. 


D. Entropy production using path probabilities: 
Fluctuation theorems 


EP along stochastic trajectories of a non-equilibrium 
system may also be estimated by using the inequality of 
probabilities of time-forward trajectories, and conjugate 
trajectories under suitably time-reversed protocol. We 
consider the time evolution of a macrospin from t = 0 
to To through a path X = [m(t),H(t)] where H(t) acts 
as control parameter, the functional form of which gives 
a specific protocol. Let us divide the path into i = 
1,2,... ,N segments, each of time-interval St such that 
NSt = tq. The transition probability pf (m', t + (5t|m, t) 
on 7-th infinitesimal segment is governed by the Gaussian 
random noise at i-th instant obeying probability distri¬ 
bution P(hi) = exp{—Sthf/4:D). Denoting 

Eq.(l2]) as rh = $(m(t),H(t)), the transition probability 
on 7-th segment pf = ^ dh.iP{h.i)5{'mi — $i), where 

the Jacobian of transformation at 7-th instant of time 
jy = det[ (9h/9m)i ]. Using Stratonovich discretiza¬ 
tion, one can show 


JA 


1 

St 


St dT{mi) 

2 dmi 


(16) 


where iF(mj) = (m x H)^ -|- 77[(a — &777^)m]j -f gili. Note 
that 9(m X H)i/dmi = 0, and d'Ri/dmi = 0. Thus the 
operative part of tF(mi) in the above relation is the ef¬ 
fective field contribution H™* = 77[(a — brrA)m]i. Ea. dTBl) 
can be expressed as 


JA = 


1 

Tt 


St 9H‘"^(m,) 

2 dmi 


(17) 


The probability of a complete trajectory is = 



Similarly, the conjugate trajectory under time-reversal 
may be discretized, and the probability of such complete 
trajectories may be expressed as V- = Y\^=iPi ■ There 
exists various possibilities to choose conjugate trajecto¬ 
ries under time-reversed protocol |29l.l42l l43l|. The conju¬ 
gate trajectory must be carefully chosen so that the ratio 
P-i_ /P_ serves as a measure of irreversibility of the pro¬ 
cess, and as a result characterizes EP in the surrounding 
environment. 

Under time-reversal, H and m changes sign simulta¬ 
neously. The corresponding conjugate trajectory is de¬ 
noted by JP! = [—m(ro — t), —H(ro — t)]. This is sim¬ 
ilar to requirement of reversal of external flow direc¬ 
tion in Ref. [d^, under time reversal. The probabil¬ 
ity of time-reversed trajectory P_ = Y\i=iPA ^ where 
p~ = J~ f dhiP(hi) S(mi - 4>i(ro - t )). It is easy to 
see from Eq. (113 that JA = ■ After some algebra, 

one obtains the ratio of the two probabilities of forward 
and reverse paths ^ = exp(Asr/^B), where 


As,- 


AL 

D 



dtn^^.m = - 


Ag 

ksT' 


(18) 
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Note that the expression of Asr presented above agrees 
with the EP given in Ea. dTK]) . Let us now assume that 
So and S£ are stochastic entropies of the system corre¬ 
sponding to its initial and final steady states respec¬ 
tively. So, the change in stochastic system entropy is 
given by As = s^ - sq = ks \n{Po/Pe) where Po(mo, Hq) 
and P^(mf, H^) are distribution functions of these micro¬ 
states. 

As we have shown above, the change in reservoir en¬ 
tropy depends on the trajectory and is given by Asr = 
ln(P+/P-). Thus the total entropy change 

Ast = ks In = As + Asr- (19) 

This immediately implies an integral fluctuation theo¬ 
rem (IFT) — I Note that in deriving 

IFT, J2x = Sx used, as the Jacobian of trans¬ 
formation from time-forward path X to time-reversed 
path is unity [4l|. Further, in a steady state, the 
total entropy change along a time-forward path As{ 
is equal and opposite to that along the time-reversed 
path, As)'(AI) = —As({X). Using this, and Ea. lfT^ 
one obtains the following detailed fluctuation theorem 
(DFT) [Ullll for probability distribution of EP p(Asi) 
as 


p(Ast) = e^*‘/'=®p(-Ast)- (20) 

Using the definition Asr = —Aq/T the IFT 

(exp(—Ast/fcs)) = 1 can be expressed as 

{exp{PAq - As/ks)) = 1. (21) 

This is equivalent to Jarzynski relation, for transforma¬ 
tions between non-equilibrium steady states Due 

to Jensen inequality, this implies T{As) > (Aq). For an 
infinitesimally slow variation of H(t), the equality holds, 
i.e., the steady state change in system entropy can be 
evaluated in terms of {Aq) ss (—m • AH). For a time- 
independent external field, one reaches an equilibrium 
steady state with (Aq) = 0, and (As) = 0. 


E. Other possibilities of conjugate trajectories 

Let us now consider, three other possibilities of choos¬ 
ing conjugate trajectories, such that one obtains entropy 
like quantities that obey DFT [i^. First, assume con¬ 
jugate trajectories such that time forward protocol of 
H(t) traces back itself under time reversal. The cor¬ 
responding conjugate trajectory is denoted by = 
[m(ro — t), H(ro — t)] where m and H do not change sign. 
The probability of such conjugate trajectories is denoted 
by Then the ratio of probabilities of time-forward 

and conjugate trajectories is P+/Pi^^ = exp(Asr^^/fcB) 


where 

Asl^^ 


ks 


D Jo dJ, 

— [ASt. J- ASgyro] ; 

Kb 


dt N • m 


( 22 ) 


where Asgyro = Awgyro/T with Awgyro = (1/d) / dlN-rh 
being the gyroscopic work done on the system due to spin 
torque. One obtains the DFT 




where As^^ = As -|- Asr^^. Numerical simulation of 
macrospins with constant amplitude m has been used in 
Ref. [i^ to obtain the probability distribution p(As(^^), 
which obeys DFT. This form of DFT may be interpreted 
as follows. One can define As* = As — Aq/T, and rewrite 
the DFT as. 


p(Ast, AWgyro) ^ + 

pi Ast, Arcgyro) 


In a steady state, ignoring As with respect to Aq/T, this 
relation leads to 


pj Aq, Awgyro) _ ^_i 3 (Aq-Aw,,^ro) ^ (’ 23 ) 

p(^9i ~AWgyro) 

This equality is closely related to the fluctuation theorem 
for heat engines , and was used in Ref. in the 

context of an isothermal engine absorbing heat Aq and 
performing work ArCgy^o via spin torque. 

The Jacobian of transformation from time forward tra¬ 
jectory X and the conjugate trajectory is unity. This 
leads to the IFT which by Jensen’s 

inequality gives (As|^^) > 0, a result equivalent to the 
second law of thermodynamics. The IFT obtained from 
Eg. (1^51) has the form = 1, which after 

Jensen’s inequality gives (Awgyro)/(Aq) < 1, meaning 
average work does not exceed average heat. Note that 
the torque N is associated with the reversible part of 
probability current J*'®'' = NP, and thus does not con¬ 
tribute to heat flux. However, it still contributes towards 
an entropy like term AsE^ that gives total entropy As^^^ 
obeying DFT and IFT. 

Next we assume that m alone changes sign along the 
conjugate trajectories so that they are described by A^ = 

[—m(To — t), H(to — t)]. We denote the path probabilities 

( 2 ) 

along such conjugate trajectories by P_ . Then the ratio 
= exp(Asr^V^s) where 


ks 


1 r° 

— / dTN.rh 

D Jo 


As 


gyro 


ks 


(24) 


Again, Asp^ = As-I-As/^^ obeys the DFT. However, the 
Jacobian of transformation from A to AI is not unity, 
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and the IFT is not obeyed by this quantity. This is ex- 
(21 

pected, as Asr depends only on N, which is associated 
with reversible probability current, and should not give 
rise to second law like inequality. 

The third alternative is to consider conjugate trajec¬ 
tories in which H alone changes sign, i.e., = [m(ro — 

t),—H(ro — t)]. Denoting the probability of conjugate 
trajectory one obtains V+j'P^"’ = 1, i.e., the corre- 

spending stochastic EP in the reservoir Asr = 0. 

The EP in reservoir associated with dissipated heat 
Asr = — Asr^\ Note that the amplitude of magne¬ 

tization can be approximated to be constant, for samples 
with Curie temperature much larger than room temper¬ 
ature. In such cases, the stochastic Langevin dynamics 
can be described as diffusion of a particle under suitable 
torque due to external field [i^. In the spherical polar 
coordinates, macrospin orientation {0, (j)) may be treated 
as even functions under time reversal. As a result one 
obtains an expression of entropy, which is equivalent to 
Asr^^ involving a gyroscopic term Asgyro. The proba¬ 
bility distribution of total EP p(AS(^^) has been shown 
to obey DFT. Of course, even within that restricted dy¬ 
namics, if one considers = [—m(ro — t), —H(ro — t)] as 
the conjugate trajectory, one obtains EP in the reservoir 
Asr = —Aq/T, as is shown in the appendix of Ref. [43| . 

Among all possible prescriptions for constructing 
stochastic trajectories, the definition of Asr in Ea. dT^ 
obtained by tracing back the time-reversed trajectory 
directly utilizing reverse protocol of H(t), such that, 
= [—m(ro — t), —H(ro — t)] leads to the expression Sr 
in Eg. lTSl) obtained from Eokker-Planck equation. The 
ratio of probabilities of time forward, and time reversed 
trajectories gives unity in presence of time reversal sym¬ 
metry. Thus any other value of this ratio gives a mea¬ 
sure of breaking of time-reversal symmetry, and thus the 
EP. Note that the derivation of Sr in Eg. dTSl) depends 
only on the dynamics, not on any particular definition 


of conjugate trajectory. Such definitions were used as 
mathematical construct to derive fluctuation theorems. 


IV. SUMMARY 

We studied stochastic thermodynamics for a macrospin 
of fluctuating amplitude and direction of magnetization 
subjected to external magnetic field. We considered 
a generalized Langevin dynamics of macrospins, taking 
into account (i) a stochastic rotational dynamics of the 
magnetization and (ii) its longitudinal fluctuations, (iii) 
a mean field approximation of the interaction between 
spins within the macrospin cluster, and an external mag¬ 
netic field. We obtained several possible fluctuation the¬ 
orems for entropy-like quantities found from using differ¬ 
ent choices of conjugate trajectories under time reversal. 
Only one of the possible choices gave As^ = —Aq/T, 
the entropy production (EP) in the reservoir due to dis¬ 
sipated heat —Aq, that agrees with the expression one 
obtains from Eokker-Planck equation. A second entropy 
like quantity As|^^ = As — Aq/T -|- Awgyro/T, where 
Awgyro is the rotational work done on the macrospin due 
to magnetic field induced spin torque, also obeys fluc¬ 
tuation theorems. The heat dissipation and gyroscopic 
work done, can be measured separately in experiments 
on macrospins, and our predictions regarding fluctuation 
theorems can be tested. 
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